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ABSTRACT 


Efficient  numerical  programs  of  great  generality  are  uti- 
lized for  steady  forced  transfer  problems  with  constant  heat  input 
per  unit  length.  Fully  developed  coupled  velocity  and  convective 
thermal  fields  for  long  straight  ducts  of  arbitrary  cross  sections 
are  obtained  by  direct  finite-difference  solution  of  the  governing 
momentum  and  energy  equations.  Arbitrary  Dirichlet  temperature 
boundary  conditions  are  considered.  Example  solutions  of  simple 
sections  are  compared  favorably  with  analytic  solutions.  Other 
examples  Include  extended  surface  sections  more  typical  of  some 
contemporary  heat  exchange  channels.  Application  to  heat  transfer 
is  discussed. 
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I . INTRODUCTION 


Prediction  of  forced  thermal  transport  in  long  straight 
ducts  of  various  cross  sections  (Fig.  1)  has  long  been  sought  theo- 
retically since  many  shapes  are  used  in  different  cooling  arrays. 
The  possibilities  of  closed  analytic  solution  have  not  yet  been 
completely  exhausted  for  fully  developed  thermal  fields,  but  the 
complexity  of  thermal  field  solution  tends  to  increase  as  the  sec- 
tion shape  deviates  from  a simple  circle  or  a planar  gap.  This  is 
due  to  the  complexity  of  the  fully  developed  laminar  velocity  field 
in  the  duct.  The  convective  heat  transfer  is  controlled  by  this 
velocity  field  so  the  resulting  thermal  field,  when  it  can  be 
found  analytically,  is  also  a complicated  solution  in  general.  In 
principle,  an  exact  analytic  solution  is  desirable,  but  in  practi- 
cal application  extensive  computer  calculations  are  often  required 
to  evaluate  quantities  that  can  be  compared  with  usual  average  mea- 
surements, much  less  point-by-point  measurements.  Many  duct  shapes 
are  intractable  to  exact  analysis.  This  report  shows  how  to  employ 
the  computer  in  direct  fashion  to  obtain  accurate  solutions  for 
arbitrarily  shaped  ducts  and  various  Dirichlet  boundary  conditions. 

A thermal  convective  field  solution  will  depend  on  the 
assigned  thermal  boundary  conditions.  Constant  heat  input  (or 
withdrawal)  per  unit  duct  length  is  hypothesized  as  it  is  compati- 
ble with  the  concept  of  a fully  developed  field  (below).  Various 
temperature  conditions  at  the  cross-section  perimeter  will  be  con- 
sidered . 


V 


Fig.  1 Long  Straight  Ducts  with  Irregular  Cross  Section 
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A fully  developed  steady-state  temperature  field  has  been 
assumed.  This  implies  a great  distance  from  the  section  where  the 
temperature  conditions  were  first  imposed  (just  as  the  fully  devel- 
oped velocity  field  implies  a long  passage  from  the  entrance).  In 
engineering  parlance,  the  theory  applies  properly  to  lengths 
greater  than  the  hydrodynamic  and  thermal  entrance  lengths.  The 
theoretical  velocity  field  is  rigorously  a condition  of  minimum 
pressure  drop  per  unit  length,  and  it  appears  that  the  correspond- 
ing thermal  field  also  represents  a minimum  Nusselt  number  situa- 
tion (see  Ref.  1,  especially  Discussion  by  Reynolds).  Together 
the  fields  provide  far  downstream  asymptotic  limits  for  developing 
fields.  Extensive  theoretical  calculations  of  development  are  be- 
yond the  range  of  the  present  report  but  are  within  today's  com- 
puter capability. 

The  equations  for  the  physical  fields  studied  here  all  in- 
volve the  well-known  Laplace  operator.  The  finite-difference 
numerical  solutions  with  constant  peripheral  temperature  are  com- 
pared to  exact  analytic  results  (some  new)  for  the  same  sections 
to  verify  their  accuracy.  The  method  is  extended  to  more  arbi- 
trary sections  with  extended  surfaces  resembling  some  contemporary 
heat  exchange  channels  for  which  no  analytic  solutions  have  yet 
been  demonstrated.  Other  temperature  boundary  conditions  also  are 
considered.  The  calculation  of  heat  transfer  quantities  used  in 
experimental  measurements  is  straightforward  and  efficient.  The 
techniques  and  conclusions  are  just  as  relevant  to  corresponding 
mass  transport  problems  (Ref.  2). 


Ref.  1.  E.  M.  Sparrow  and  R.  Siegel,  "A  Variational  Method 
for  Fully  Developed  Laminar  Heat  Transfer  in  Ducts,"  Trans.  ASME. 
Vol.  79C,  1959,  pp.  157-167. 

Ref.  2.  W.  M.  Kays,  Convective  Heat  and  Mass  Transfer. 
McGraw  Hill,  New  York,  1966. 
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II  EQUATIONS  AND  BOUNDARY  CONDITIONS 


The  physical  problems  will  be  idealized  to  constant  prop- 
erty fluid.  The  fluid  is  incompressible  and  has  Newtonian  vis- 
cosity y and  density  p.  Frictional  heating  will  be  neglected. 
The  thermal  properties,  diffusivity  a,  conductivity  k,  and  spe- 
cific heat  Cp  of  the  fluid  are  constant.  The  geometries  will  be 
straight  and  infinitely  long  in  the  z direction  (Fig.  1).  The 
steady  laminar  flow  is  fully  developed,  i.e.,  there  is  only  one 
(z-ward)  velocity  component  u,  and  the  longitudinal  pressure 
gradient  is  constant  and  uniform  over  the  section.  The  single 
component  momentum  equation  reduces  to 


2.  i*£ 

y dz 


(1) 


Normalizing  by  the  duct  half-height  in  the  central  plane,  b, 
Eq.  (1)  may  be  rewritten  in  nondimensional  form 


2 * 

V u - -1 

where 

* u 
u = — 

kl  i£ 

u dz 


(la) 


(lb) 


Many  solutions  to  this  equation  (with  the  nonslip  viscous 
boundary  condition  u*  * 0 on  the  duct  wall)  for  various  cross- 
sectional  shapes  are  available,  some  analytical,  some  numerical. 

In  principle,  the  velocity  field  for  any  shaped  duct  can  be  ob- 
tained via  transformation  of  the  cross  section  to  a circle  (Ref.  3). 
In  practice,  it  is  easier  to  work  in  a Cartesian  coordinate  system 
and  provide  the  solution  numerically  (Ref.  4).  The  finite-differ- 
ence systems  are  all  solved  directly. 


Ref.  3.  L.  N.  Tao,  "On  Some  Laminar  Forced-Convection  Prob- 
lems," Trans.  ASME.  Vol.  83C,  1961,  pp.  446-472. 

Ref.  4.  V.  O'Brien,  "Steady  and  Unsteady  Flow  in  Non-Circu- 
lar Straight  Ducts"  (accepted  by  ASME  J.  Appl.  Mech. , 1976). 
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A fully  developed  temperature  field  t with  3t/3z  constant 
with  respect  to  both  flow  length  and  cross-sectional  position  has 
been  specified  earlier  for  constant  heat  input  per  unit  length  and 
constant  wall  temperature.  We  note  these  are  idealizations,  and  it 
is  permissible  to  extend  the  concept  of  a fully  developed  field  to 
the  cases  where  the  heat  input  or  temperature  vary  around  the 
perimeter.  Using  the  subscripts  w for  wall  and  ref  for  a refer- 
ence temperature  at  the  section,  along  with  an  overbar  for  sec- 
tional average,  a generalized  fully  developed  thermal  field  re- 
quires 


Often  the  reference  temperature  is  (tw  - t^)  where  the  bulk  (mixed 
mean  or  "mixing-cup")  temperature  for  the  cross-sectional  area  A 

is  defined:  t^  = JAut  dxdy.  It  is  assumed  that  tw  and  tb  are 

Au 

functions  of  z only.  Carrying  out  the  differentiation, 


Since  for  the  case  of  constant  heat  input  neither  of  these  right- 
hand  terms  is  equal  to  zero,  8t/3z  must  equal  dtw/dz  in  order  to 
cancel  the  dependence  on  cross-sectional  position.  Then  dtw/dz 
« dtb/dz,  and  by  doing  an  energy  balance  per  unit  length  of  chan- 
nel dtb/dz  is  a constant.  Thus 


LL 

3 z 


dt  dt, 

w b 

dz  dz 


Constant 


2_ 

(A)  pc  u 
P 


(4) 


where  Q « volume  heat  addition  per  unit  channel  length.  If  tw  is 
absolutely  constant,  on  the  other  hand,  and  not  dependent  on  z. 


tlf  the  quantity  is  defined  for  every  point  of  the  section,  the 
overbar  represents  an  area  average;  if  the  quantity  is  defined 
on  the  wall,  the  overbar  represents  a perimeter  (circumferential) 
average. 
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t - t dt,  n 

— = — ~ x — — « Constant  = ^ — 3 (5) 

3z  Cw  " Cb  dZ  (A)pcp  u 

Then  dtb/dz  is  not  constant,  and  a numerical  solution  is  more  in- 
volved (see  below). 

It  is  customary  in  forced  convection  problems  to  neglect 
the  heat  diffusion  in  the  flow  direction  and  natural  (free)  con- 
vection due  to  gravity  force.  We  also  assume  no  internal  heat 
sources  in  the  fluid.  The  steady  thermal  transport  equation  for 
a fully  developed  temperature  field  is  a form  of  the  energy  equa- 
tion (Ref.  1). 


_u  3t_ 
a 3z 


(6) 


where  a = k/pcp  is  the  thermal  diffusivity.  Because  3t/3z  is  con- 
stant, a normalized  equation  for  the  thermal  field  appropriate  for 
the  constant  heat  input  per  unit  length  can  be  written 

2**2*—  ,,  . 

7 t - u = V (t  - t ) , (6a) 

w 

where  t*  has  been  nondimensionalized , via  Eqs.  (lb)  and  (4),  as 


-_E (6b) 

b dp  3t 
pa  dz  3z 


Successful  development  of  practical  heat  transfer  theory 
requires  boundary  conditions  appropriate  to  the  physical  situations. 
Most  prior  efforts  have  involved  homogeneous  boundary  conditions, 
which  has  simplified  the  mathematical  algebraic  analysis.  However, 
taking  into  account  that  there  can  occur  problems  where  uniform 
boundary  values  are  not  valid  assumptions,  the  present  numerical 
programs  allow  arbitrary  peripheral  temperature  distributions  to 
match  a variety  of  cooling  (heating)  problems. 
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For  example,  many  applications  of  cooling  ducts  involve  an 
array  of  identical  ducts  through  a hot  material  matrix.  Unless 
they  are  lined  with  a material  of  very  high  heat  conductivity,  it 
is  unlikely  that  the  embedded  ducts  will  have  uniform  wall  tempera- 
ture. The  peripheral  temperature  will  depend  on  the  arrangement 
of  the  ducts  within  the  material  to  be  cooled.  Figure  2 shows  two 
such  arrangements.  If  the  homogeneous  matrix  contains  uniformly 
distributed  heat  sources  (e.g.,  from  radioactivity  or  exothermic 
chemical  reaction),  probable  temperature  distributions  on  the 
walls  of  the  ducts  are  indicated. 

Uniform  flux  per  unit  surface  area  at  the  wall  has  also 
been  a theoretical  assumption  that  is  unlikely  to  be  found  in  non- 
circular open  duct  forced  transfer  application.  For  extended  sur- 
faces, the  probabilities  of  uniform  temperature  or  flux  are  even 
less  (Ref.  5).  The  present  program  allows  any  assigned  tempera- 
ture distribution  to  be  applied  as  the  duct  wall  boundary  condi- 
tion (B.C.).  The  considerable  effects  of  nonuniform  temperature 
are  illustrated.  (In  principle,  one  can  solve  the  nonhomogeneous 
boundarv  value  problem  algebraically  as  well  [see  the  appendix]  , 
but  the  effects  seem  not  to  have  been  shown  even  for  relatively 
simple  shapes.) 


(a)  Rectangular  array 


(b)  Staggered  array 


Fig.  2 Different  Duct  Arrays  with  Indicated  Possible  Wall  Temperature 
Distributions.  Cross  ( + ) marks  local  peak  temperature. 


Ref.  5.  K.  A.  Gardner,  "Efficiency  of  Extended  Surface," 
Trans.  ASME,  Vol.  67,  1945,  pp.  621-627. 
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For  convenience,  homogeneous  and  nonhomogeneous  B.C.'s 
are  discussed  in  separate  sections. 


A.  ASSIGNED  WALL  TEMPERATURES 

1.  Uniform  Wall  Temperature,  t 

For  sectionally  uniform  but  linearly  varying  tw(z),  the 
homogeneous  boundary  condition  for  F.q.  (6a)  is  6*  = t*  - tw*  » 0 
on  the  duct  wall.  An  analytic  solution  for  a square  duct  based  on 
a recent  solution  of  the  velocity  field  is  given  in  the  appendix. 

It  will  be  used  to  assess  the  absolute  accuracy  of  the  present 
direct  finite-difference  solution  of  the  same  problem  (and  to  pro- 
vide an  analytic  basis  for  comparison  of  the  effect  of  other  bound- 
ary conditions).  Other  rectangular  duct  shapes  are  considered,  as 
are  curved  walls. 

Advantage  has  been  taken  of  assumed  bilateral  symmetry  to 
consider  just  one-quarter  of  the  field.  This  saves  computer  time 
but  is  not  a fundamental  requirement.  Similar  calculations  are 
possible  with  more  irregular  duct  shapes. 

2.  Variable  Wall  Temperature  t (X  , Y ) 

c w — w — w 

If  the  material  of  the  duct  walls  is  homogeneous  and 
source-free,  and  has  high  heat  conductivity,  uniform  tw*  is  proba- 
bly a valid  assumption.  For  other  cases,  temperature  distributions 
on  the  duct  section  perimeter  provide  more  general  "arbitrary" 
boundary  conditions  on  the  convective  thermal  field.  The  only 
physical  constraint  on  the  generality  is  that  the  profile  be  con- 
tinuous, particularly  at  corners  of  rectangular  ducts. 

An  analytic  solution  for  simple  temperature  profiles  on 
the  edges  of  the  square  duct  has  been  derived  for  comparison  to 
the  numerical  solution  (see  the  appendix).  Due  to  the  linearity 
of  Eq.  (6),  the  solution  consists  of  the  sum  of  the  solution  dis- 
cussed in  the  previous  section  and  a solution  satisfying  the  non- 
constant temperature  condition  on  the  perimeter. 


B.  OTHER  THEORETICAL  BOUNDARY  CONDITIONS  (PREVIOUSLY  USED  FOR 
CERTAIN  PROBLEMS) 

The  condition  that  the  duct  wall  temperature  be  an  abso- 
lute constant  independent  of  x,  y,  or  z has  also  been  used  theo- 
retically for  simple  duct  problems.  This  is  a very  special  case 
For  tw  constant,  Eq.  (5)  substituted  into  Eq.  (6)  reduces  it  to 
(approximately) 
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(V  + b'u  )t  * a'u 


(a' ,b'  are  constants). 


(7) 


This  is  a more  complex  equation  than  Eq.  (6a)  and,  although  the 
operator  (V^  + b'u*)  can  be  iterated  numerically  (Ref.  6),  such 
solutions  for  noncircular  ducts  will  not  be  considered  here.  The 
wall  boundary  condition  is,  of  course,  t - tw.  Previously  pub- 
lished results  for  simple  shapes  show  this  boundary  condition  pro- 
vides less  efficient  average  heat  transfer  (i.e.,  lower  Nusselt 
numbers)  than  variable  tw(z)  (Ref.  2,  Table  8-2). 

Some  fully  developed  temperature  fields  have  been  analyzed 
with  assigned  piecewise  uniform  flux  boundary  conditions  (Refs.  1 
and  7).  Rectangular  duct  solutions  were  obtained  by  separating  the 
total  temperature  solution  into  two  portions.  The  first  takes  an 
analytic  form  exactly  satisfying  the  assigned  flux  condition  (gradi- 
ent or  Neumann  condition)  and  Eq.  (6)  — call  it  tf.  This  reduced 
the  second  portion  of  the  solution  — call  it  to  - to  a zero  flux  con- 
dition on  the  boundary  which  had  to  satisfy  (in  the  present  notation) 


* 


* — * 
u - u . 


(8) 


The  rectangular  solution  tg*  Is  independent  of  the  actual  flux  con- 
ditions (such  as  given  relative  flux  on  alternate  sides).  For  less 
regular  sections  (nonrectangular , noncircular),  it  is  more  diffi- 
cult to  find  the  exact  (analytic)  flux  portion,  tf,  for  uniform 
flux  (or  other  assigned  distributions).  This  portion  is  constrained 
(mathematically  and  physically)  to  satisfy 

2 * — * 

V tf  - u (9) 


Ref.  6.  S.  H.  Clark  and  W.  M.  Kays,  "Laminar-Flow  Forced 
Convection  in  Rectangular  Tubes,"  Trans.  ASME,  Vol.  75,  1953,  pp. 
859-866. 

Ref.  7.  J.  M.  Savino  and  R.  Siegel,  "Laminar  Forced  Convec- 
tion in  Rectangular  Channels  with  Unequal  Heat  Addition  on  Adjacent 
Sides,"  Int.  J.  Heat  Mass  Transf,.  Vol.  7,  1964,  pp.  733-741. 
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Table  14.02  In  Ref.  8 may  be  useful  for  separable  coordinate  sys- 
tems, but  the  flux  boundary  value  problem  is  more  complicated 
analytically  and  numerically  than  for  Dirichlet  temperature  con- 
ditions. With  uniform  flux  boundary  conditions  on  noncircular 
ducts,  the  wall  temperature  varies  around  the  perimeter  (Refs.  1 
and  7). 

The  so-called  radiation  boundary  condition  3t/3n  ■ Xt 
(n  ■ normal  to  surface)  has  also  been  used  as  a theoretical  con- 
vective field  boundary  condition  for  one-dimensional  flows.  Ex- 
tending this  to  noncircular  ducts  implies  3/3n(fen  t)  ■ X or 
t • tw  e^n  uniformly  over  the  fluid  region  enclosing  the  duct 
wall.  Generally  3u*/3n  varies  around  the  perimeter  while  at  the 
wall  3/3n(V2t*)  - 33t*/3n3  ■ constant  by  the  hypothesized  B.C. 

So  this  condition  does  not  seem  appropriate  for  fully  developed 
fields  in  irregular  duct  shapes. 
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III.  NUMERICAL  SOLUTIONS 


This  section  presents  a survey  of  the  kinds  of  steady  de- 
veloped convective  field  solutions  that  can  be  obtained  directly 
with  computer  programs  using  the  Cartesian  finite-difference  formu- 
lation of  the  Laplace  operator,  V^.  The  numerical  schemes  them- 
selves are  variations  of  Ref.  9 and  were  discussed  very  briefly  in 
the  appendix  of  Ref.  4 by  Dr.  Louis  W.  Ehrlich  who  programmed  the 
various  problem  sets.  The  fully  developed  velocity  field  is  ob- 
viously an  integral  part  of  each  solution,  but  these  problems  have 
been  considered  at  length  in  Ref.  4.  Without  repeating  the  dis- 
cussion, the  main  point  is  that  the  numerical  velocity  solutions 
are  quite  accurate^  As  pointed  out  above,  the  thermal  field  solu- 
tions will  depend  upon  the  thermal  boundary  conditions  assigned  on 
the  duct  walls. 


A.  UNIFORM  PERIPHERAL  WALL  TEMPERATURE 
1.  Rectangles 

The  complete  numerical  solution  with  h = 1/32  for  the  nor- 
malized temperature  difference  6*  - t*  - tw*  with  6*  ' 0 on  the 
duct  walls  for  a square  is  illustrated  in  Table  1 by  the  center- 
plane  temperature  distributions.  Both  the  absolute  temperature 
and  that  relative  to  the  calculated  central  temperature  0C*  are 
given  for  comparison  to  the  analytic  solution.  The  agreement  is 
very  good,  and  if  higher  digit  accuracy  is  desired  a finer  mesh 
solution  is  easily  attainable  at  a slightly  higher  cost  in  computer 
time.  The  cost  of  obtaining  the  velocity  and  temperature  values  at 
nearly  a thousand  interior  points,  h-^  = (32) was  $2.13  on  the 
IBM  360/91  computer.  The  cost  of  evaluating  the  analytic  solutions 
for  16  points  on  the  central  plane  was  greater  on  the  same  com- 
puter. 


Values  of  the  temperature  distribution  in  the  narrow  cen- 
tral plane  of  other  rectangular  ducts  are  also  given  in  Table  1. 
Some  comparison  values  from  a previous  pioneering  numerical  solu- 
tion (Ref.  6)  are  also  shown.  It  is  evident  that  the  profile 


Ref.  9.  L.  W.  Ehrlich,  "Solving  the  Biharmonic  Equation  as 
Coupled  Finite  Difference  Equations,"  SIAM  J.  Num.  Anal..  Vol.  8, 
1971,  pp.  278-287. 

tlllustrative  examples  showed  errors  in  the  fourth  or  fifth  signifi- 
cant digit. 
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Table  1 

Centerplane  Temperature  Profiles  in  Rectangular  Ducts 


i 

Square  (a/b  - 1) 

Planar  gap  (a/b  ■ °°) 

1 

Absolute 

A A 

0 (0,Y)  or  6 (X,0) 

Relative 

(6c*)-1  6* 

Relative 
(Sc*)"1  9* 

Y,(X) 

Numerical 

Analytic 

Numerical 

Analytic 

Analytic 

0 

-0.064990 

-0.064998 

1.0 

1.0 

1.0 

0.125 

-0.063841 

-0.063848 

0.9823 

0.9823 

0.9813 

0.25 

-0.060413 

-0.060418 

0.9296 

0.9296 

0.9258 

0.375 

-0.054766 

-0.054769 

0.8427 

0.8426 

0.8352 

0.5 

-0.047010 

-0.047011 

0.7233 

0.7233 

0.7125 

0.625 

-0.037324 

-0.037322 

0.5743 

0.5742 

0.5618 

0.75 

-0.025976 

-0.025972 

0 3997 

0.3996 

0.3883 

0.875 

-0.013353 

-0.013348 

0.2055 

0.2054 

0.1985 

1.0 

0 

0 

0 

0 

0 

Rectangle  (a/b  - 2) 


Narrow  Profile 

Wide  Profile 

Absolute 

Relative 

X 

Absolute 

Relative 

-0.162114 

-0.162114 

1.0 

-0.1626+ 

1 • u 

U 

-0.1626+ 

-0.159106 

0.9815 

0.25 

-0.159926 

0.9865 

-0.150172 

0.9263 

0.5 

-0.153244 

0.9453 

-0.135575 

0.8363 

0.75 

-0.141733 

0.8743 

-0.115761 

0.7141 

1.0 

-0.124893 

-0.1256+ 

0.7704 

-0.091364 

0.5636 

1.25 

-0.102199 

0.6304 

-0.063212 

0.3899 

1.5 

-0.073340 

0.4524 

-0.032340 

0.1995 

1.75 

-0.038678 

0.2386 

0 

0 

2 

0 

0 

tRenormalized  value  from  Clark  and  Kays  relaxation  solu- 
tion (Ref.  6). 
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shapes  are  very  nearly  identical  in  the  narrow  central  plane;  the 
major  variation  with  aspect  ratio  a/b  is  the  central  value,  0C*. 

Temperature  distributions  on  the  wider  central  plane  are 
illustrated  by  Fig.  3.  The  shape  of  the  profile  becomes  flatter 
as  the  width-to-height  ratio  ■ a/b  approaches  infinity  (Fig.  3a). 
If  the  relative  results  are  plotted  in  "stretched  coordinates," 

6 /0C*  versus  X'  * (b/a)X,  it  is  more  evident  that  the  curves  are 
not  similar  but  vary  with  a/b  (Fig.  3b).  The  variation  along  the 
diagonals  of  the  rectangles  is  also  illustrated  in  Fig.  4,  where 
S'  *=  2_is[Y2  + (X')2]^  is  another  stretched  coordinate.  Also  in- 
cluded in  Fig.  4 is  the  corresponding  profile  from  a square  duct 
approximate  thermal  solution  via  a variational  approach  (Ref.  1). 

2 . Ellipses  and  Other  Curved  Shapes 

Similar  numerical  results  are  given  for  circular  and 
ellipsoidal  ducts  in  Table  2 and  Fig.  5.  The  finite-difference 
calculations  were  carried  out  in  Cartesian  coordinates.  The  com- 
parison with  the  conformal  analytic  results  (Ref.  3)  is  quite  rea- 


(b)  Relative  profile  versus 
"stretched  coordinate,"  X' 


Fig.  3 Wide  Centerplane  Temperature  Profiles  in  Rectangular  Ducts  with 
Uniform  Peripheral  Temperature  for  Several  a/b 
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Fig.  4 Absolute  Diagonal  Temperature  Profiles  for  Rectangular  Ducts  with 
Uniform  Peripheral  Temperature  as  a Function  of  a/b 
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t v ± X' 


(a)  Relative  narrow  profiles  (b)  Relative  wide  profiles 

Fig.  5 Centerplane  Temperature  Profiles  in  Elliptical  Ducts  as  a Function  of  a/b 


sonable.  (On  the  scale  of  Fig.  5,  the  numerical  and  analytic  re- 
sults cannot  be  distinguished.)  The  narrow  central-plane  relative 
profiles,  6*/ec*,  are  essentially  the  same,  independent  of  the 
semimajor  axis  ratio  (a/b),  and  very  similar  to  those  for  the 
rectangles  shown  on  Table  1.  However,  as  a/b  ® the  wide  center- 
plane  profiles  tend  to  bell-shapes  more  like  the  diagonal  profiles 
in  the  rectangles. 

The  "short-legged"  formulas  applied  to  Cartesian  coordi- 
nates to  get  these  results  for  the  conic  sections  can  be  applied 
just  as  well  to  other  curved  contours.  Any  combination  of  flat 
and  circularly  filleted  arcs  could  be  handled  just  as  readily.  It 
is  not  necessary  to  have  an  analytic  description  of  the  wall  con- 
tours. The  well  coordinates  can  be  read  in  from  rectangular  grid 
measurements. 
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3.  Extended  Surfaces 

Heat  exchange  ducts  often  use  fins  or  other  extensions  to 
increase  the  relative  surface  (Refs.  5 and  10).  The  numerical 
scheme  can  be  used  for  contours  including  extended  surfaces  as 
illustrated  in  Fig.  6 for  uniform  peripheral  temperature  on  a 
square  section  with  rectangular  fins.  As  might  be  expected,  the 
isotherms  tend  to  take  on  the  shape  of  the  wall  contours,  but  they 
are  not  identical  with  isovels  of  the  flow  field.  The  peak  veloc- 
ity (marked  by  a cross)  has  been  displaced  from  the  centerline  of 
the  duct,  but  the  greatest  temperature  difference,  peak  0*,  still 
occurs  there. 

The  effect  of  fin  length  is  illustrated  for  narrower  fins 
in  the  same  square  duct  by  Fig.  7.  Here  the  emphasis  has  been 
placed  on  the  temperature  gradient  along  the  perimeter  of  the  duct 
section.  (The  no-fin  case  is  included  for  contrast.)  Again,  the 
peak  velocity  in  the  finned  duct  is  indicated  for  each  case,  and 
it  should  be  noted  that,  as  the  fin  length  increases,  the  peak 
velocity  moves  away  from  the  duct  centerline  and  is  replaced  by 
twin  peaks  that  migrate  toward  the  centers  of  two  almost  rectangu- 
lar ducts.  (In  the  limit  when  the  fins  touch,  there  are  two  sepa- 
rate rectangular  ducts.) 


B.  ASSIGNED  WALL  TEMPERATURE  DISTRIBUTION 

1.  Rectangles  ^ 

The  square  duct  solution  with  assigned  parabolic  tempera- 
ture distributions  along  the  walls  (Fig.  8)  was  obtained  by  the 
direct  numerical  technique.  The  average  0W*  « 0.  The  results  are 
dependent  on  a nondimensional  parameter  e (see  the  appendix)  that 
characterizes  the  magnitude  of  the  nonuniformity  but  can  be  posi- 
tive or  negative.  The  difference  from  the  uniform  0W*  = 0 problem 
can  be  considerable  if  |e|  is  not  small.  Isotherms  are  compared 

in  Fig.  8 and  (normal)  wall  temperature  gradients  (or  flux)  in  1 

Fig.  9.  (The  accuracy  of  the  finite-difference  calculation  I 

(h  ■ 1/16)  with  asslgn^d^wall  distribution  is  verified  by  the  cen- 
ter value  0C*  - tc  - tw  ; analytically  8C*  « -0.5288,  numerically  _ 

0C  * -0.5315,  about  0.5%  deviation.)  t 


• 10. M. H. Hu  and  Y.  P.  Chang,  "Optimization  of  Finned 
Tubes  for  Heat  Transfer  in  Laminar  Flow,"  Trans.  ASME.  Vol.  95C. 
1973,  pp.  332-338. 
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0 

(a)  Isovels 


V 


(b)  Isotherms 


Fig.  6 Velocity  and  Temperature  Fields  for  a Square  Duct  with  Extended 
Surfaces  (upper  half  shown) 


Arbitrary  scale 
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C Fin  B A pjn  AD  E 

tip  side  Rest  of  wall 


Fig.  7 Normal  Derivatives,  d0*/3n,  along  the  Duct  Surface  with  Uniform 
Peripheral  Temperature  for  Three  Narrow  Finned  Configurations 
(one  quadrant  shown) 


I 
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y 


(a)  Isotherms 


Fig.  8 Thermal  Field  for  a Square  Duct  with  Parabolic  Temperature  Distributions  on 
the  Perimeter  (peaks  at  corners),  (a)  and  (b)  Wall  temperature  distributions; 
(c)  isotherms. 
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2.  Extended  Surfaces 


The  numerical  solutions  with  assigned  wall  values  can  be 
extended  to  curved  contours  and  extended  surfaces.  Fins  with  tem- 
perature distributions  (Fig.  10a)  combined  with  duct  walls  with  a 
uniform  value  give  distorted  thermal  fields  and  very  different 
flux  values  along  the  surface  (Fig.  10b)  when  compared  to  totally 
uniform  peripheral  temperatures.  The  difference  in  magnitude  of 
the  normal  derivative  is  so  large  the  uniform  case  appears  as  es- 
sentially zero.  The  parameter  6,  which  is  the  value  at  the  tip 
comer,  characterizes  the  effect  for  these  similar  temperature  dis- 
tributions. 
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Side  of  fin  End  of  fin 

(a)  Wall  temperature  distributions  on  the  fins 


Fig.  10  A Finned  Square  Duct  with  Various  Similar  Symmetric  Temperature 
Distributions  on  the  Fins 
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IV.  HEAT  TRANSFER  APPLICATIONS 


A.  DEFINITIONS 
1.  Average  Heat  Transfer 

Average  thermal  fields  are  usually  defined  for  engineering 
use  in  terms  of  a "mixed  mean"  or  bulk  temperature  tfc,  instead  of 
a simple  mean  temperature.  An  average  heat  transfer  conductance  h 
a a often  based  on  the  total  heat  transfer  rate  to  the  fluid  and  the 
temperature  difference  between  the  wall  of  the  duct  and  the  bulk 
temperature  (Refs.  1 through  3). 


h 


S 

(heated  perimeter)  • (t^  - tw> 


This  is  usually  expressed  in  terms  of  a nondimensional  (average) 
Nusselt  number,  Nu: 


h D 


. ~ . ,,  ,,  „ „ _ 4 cross-sectional  area 

where  De  is  the  hydraulic  diameter;  De  = — 

J » e perimeter 

■ — . In  the  present  notation,  for  open  rectangular  ducts  of  as- 
pect ratio  (a/b) , 


Nu  - 


(a/b)2 
(1  + a/b)2 


which  Includes  a Dg  of  4b  a/b)  t'iat  e<lual8  the  "geometric  diame- 

ter," 2b,  only  for  a square  duct.  This  assumes  also  that  the 
heated  perimeter  and  wetted  perimeter  are  the  same. 

Alternately  the  Nusselt  number  can  be  defined  in  terms  of 
the  geometric  diameter,  d ■ 2b  (a/b  3e  1), 
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— A 


Nu 


• = M _ IZbih  . 2u (a/b)  . rectanKle8 

" k k * (1  + a/b)  IOr  rectan«ie8 


where  average  conductance 


h = k 


Nu' 

2b 


The  length  k/h  ■ 2(Nu')~^(b)  has  been  Interpreted  (Ref.  11)  as  a 
thermal  barrier  that  "represents  the  thickness  of  a stationary 
fluid  film  having  a thermal  resistance  corresponding  to  h."  For 
fully  developed  flow  characterized  by  one  coordinate  (1-D  flow) 
and  uniform  peripheral  temperature,  namely  the  classic  circular 
duct  and  planar  gap  problems,  the  values  of  h/k  are  2.18  b“l  and 
2.06  b”l,  respectively  (see  Table  2 in  Ref.  11).  Thus  the  1-D 
thermal  barrier  layer  seems  to  be  about  one-half  the  distance  to 
the  duct  centerline. 

2.  Local  Heat  Transfer 


For  some  applications,  the  interest  is  on  the  local  sur- 
face heat-transfer  conductance  that  varies  around  the  perimeter 
of  the  noncircular  ducts  (2-D  flow).  The  surface  conductance  is 
defined  (Ref.  6)  as 


h 


loc 


where  normal  N is  X,  Y respectively  for  alternate  sides  of  an 
open  rectangular  duct.  If  is  averaged  over  the  total  perime- 
ter P,  the  mean  surface  conductance  is  ^ J h^oc(s)ds. 

oc  p oc 


Ref.  11.  R.  N.  Norris  and  D.  D.  Streid,  "Laminar-flow  Heat 
Transfer  Coefficients  for  Ducts,"  Trans.  AS ME , Vol.  62,  1940,  pp. 
525-541. 
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B.  UNIFORM  PERIPHERAL  TEMPERATURE 

1.  Heat  Transfer  in  Rectangular  Ducts 

A compilation  of  Nusselt  number  evaluations,  numerical  and 
analytic,  for  various  rectangular  ducts  is  given  in  Table  3.  Pre- 
vious numerical  results  (see  Ref.  2)  are  enclosed  in  brackets.  The 
analytic  values  have  been  reevaluated  with  series  truncated  at  a 
very  large  number  of  terms  so  there  are  some  slight  differences 
from  Ref.  12.  The  conventional  Nusselt  number,  Nu,  is  monotonic 
(increasing)  with  a/b  increasing.  However,  Nu'  has  a minimum  near 
a/b  = 2.^  Even  so  the  average  thermal  barrier  layer  varies  only 
±10%  from  a mean  value  about  (1.7)~^b  for  a/b  £ 8 (Fig.  11). 

The  variation  of  normalized  surface  conductance  ^loc/^loc 
is  illustrated  in  Fig.  12  for  the  square  duct  (compare  Fig.  2 in 
Ref.  12).  It  has  been  stated  that  a numerical  solution  is  not  suf- 
ficiently accurate  to  describe  the  variation  adequately  (Ref.  6). 
This  statement  is  usually  not  accompanied  by  a direct  comparison 
of  a fine  mesh  calculation  and  the  corresponding  analytical  eval- 
uation, as  here.  The  figure  establishes  the  engineering  utility 
of  direct  numerical  solutions  for  local  conductance  (less  than  2% 
deviation).  The  variational  solution  (Ref.  1)  is  less  accurate. 

For  further  comparison  and  physical  insight,  the  shear  profile  of 
the  velocity  field  is  also  shown  (dashed  curve).  The  variation  of 
shear  and  surface  conductance  with  distance  along  the  duct  perime- 
ter is  rather  similar.  The  peak  values  are  at  the  center,  the 
average  values  about  60%  of  the  way  to  the  corner,  and  both  curves 
fall  to  zero  at  the  corner.  Thus  the  local  thermal  barrier  layer 
(inversely  proportional  to  hioc)  increases  without  limit  as  the 
sharp  concave  corner  is  approached  and  the  velocity  shear  ap- 
proaches zero. 

2 . Heat  Transfer  in  Elliptical  Ducts 

A discussion  of  heat-transfer  results  for  elliptical  ducts 
is  found  in  Ref.  3.  There  is  no  need  to  repeat  it  here.  Rectangu- 
lar and  elliptical  sections  have  the  same  limiting  planar  gap  cen- 
terplane  temperature  distribution  as  a/b  -*•  °°,  but  the  Nu  variation 
with  a/b  is  different  (see  Fig.  4,  Ref.  5).  One-half  Nu'  and  its 
inverse,  the  average  thermal  barrier  layer,  for  elliptical  ducts 
rapidly  approach  limiting  values  as  a/b  increases  (Fig.  11). 


Ref.  12.  S.  M.  Marco  and  L.  S.  Han,  "A  Note  on  Limiting 
Laminar  Nusselt  Number  in  Ducts  with  Constant  Temperature  Gradient 
by  Analogy  to  Thin-Plate  Theory,"  Trans.  ASME,  Vol.  77,  1955,  pp. 
625-630. 

tMinimum  Nu'/2  ■ 1.54613  occurs  at  a/b  ■ 1.97;  maximum  barrier 
- 0.64686  b. 
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Table  3 

Nusselt  Numbers  for  Rectangular  Ducts  with 
Uniform  Peripheral  Temperature 


Aspect  ratio 
a/b 

Nu 

Nu'/2  ~ (barrier  layer)  * 

Numerical 

Analytic1' 

Numerical 

Analytic 

1.0 

3.600 

3.6079 

1.800 

1.804 

[3.63] 

(1.82) 

1.4 

[3.78] 

3.734 

(1.62) 

1.600 

1.5 

3.787 

3.790 

1.578 

1.579 

2.0 

4.131 

4.123 

1.55 

1.546 

[4.11] 

(1.54) 

3.0 

4.795 

1.598 

4.0 

5.365 

5.331 

1.68 

1.666 

[5.35] 

(1.67) 

8.0 

[6.60] 

6.498 

(1.85) 

1.828 

oo 

8.235 

8.235 

2.059 

2.059 

Note:  Bracketed  values  are  from  Kays,  Convective  Heat  and  Mass 

Transfer  (Ref.  2),  Table  8-2,  p.  117. 

tSee  the  appendix  for  the  evaluation  scheme;  there  are  slight 
differences  from  the  original  values  of  Ref.  12. 
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Fig.  1 1 Average  Thermal  Barrier  Inverse  (Nu'/2)  as  a Function  of  Aspect 
Ratio  a/b  for  Open  Rectangular  Ducts  and  Open  Ellipsoidal  Ducts 
with  Uniform  Peripheral  Temperature 
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For  elliptical  ducts,  the  shear  and  normal  temperature 
gradient  distributions  are  monotonic  on  the  quarter-perimeter. 

The  maxima  occur  at  the  narrow  centerplane  intersections , the 
minima  at  the  wide  centerplane  intersections,  the  difference  be- 
tween extrema  increasing  with  eccentricity  of  the  ellipse. 

3.  Heat  Transfer  in  Ducts  with  Extended  Surfaces 

The  finned  ducts  have  increased  perimeters  and  decreased 
areas  compared  to  the  open  square  ducts.  This  means  the  hydraulic 
diameter  De  is  always  smaller  for  the  finned  ducts.  One  gets  the 
impression  from  the  literature  that  this  decrease  in  De  is  more 
than  compensated  for  by  more  effective  heat  transfer. 

The  ordinary  Nusselt  number  for  any  noncircular  duct  is, 
in  the  present  notation. 


Nu 


2 -* 
(2A) 

Pb  S* 


The  corresponding  form  of  Nu'  is 


Nu’ 


(— ) 
w 


Values  of  P,  A,  Nu,  and  Nu'/2  are  given  in  Table  4 for  the  finned 
square  ducts  with  uniform  peripheral  temperature.  For  the  thick 
50%  fins,  the  ordinary  Nu  is  considerably  smaller,  while  the 
"thermal  barrier"  is  essentially  unchanged  from  the  open  square 
duct.  For  the  narrow  fins  of  increasing  length,  a small  distur- 
bance of  the  flat  wall  substantially  decreases  Nu,  but  as  the  fin 
gets  longer  Nu  recovers  and  exceeds  the  value  for  the  open  square 
duct.  On  the  other  hand,  Nu*  shows  less  of  a decrement  for  small 
fins  and  considerable  improvement  (Nu'  increasing,  barrier  layer 
decreasing)  as  the  fin  gets  longer.  So  by  both  criteria  only  the 
long  thin  fin  does  achieve  a more  effective  average  heat  transfer. 

The  level  of  the  temperature  gradient  along  the  unmodified 
flat  walls  is  always  decreased  by  the  presence  of  the  fins,  but 
the  fins  themselves  contribute  heavily  to  the  flux  along  the  finned 
wall.  The  distribution  of  is  similar  to  the  velocity  gradient 

distribution  (Fig.  13).  There  are  discontinuities  at  the  sharp 
corner  of  the  fin  due  to  lack  of  mathematical  normal  at  the  point. 
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Fig.  13  Variation  of  Valocity  Gradients  (shear)  and  Temperature  Gradients 

(conductance)  along  the  Perimeter  for  a Finned  Square  Duct  with  Uniform 
Peripheral  Temperature 
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Table  4 

Average  Heat-Transfer  Nusselt  Number  for 
Finned  Square  Ducts 


d/b 

Pb  1 

Ab-2 

— * 

u 

* 

«b 

Nu 

Nu’/(2) 

Wide  Fin 

0.50 

10 

3 

0.0353 

0.00594 

2.14 

1.78 

Narrow 

Fins 

0.25 

9 

3.75 

0.0927 

0.0228 

2.82 

1.70 

0.50 

10 

3.5 

0.0538 

0.00861 

3.06 

2.18 

0.75 

11 

3.25 

0.0314 

0.00285 

3.84 

3.25 

No  Fin 

0 

8 

4 

0.1404  | 

0.0390 

3.60 

1.80 

but  of  course  the  real  surfaces  are  smooth  so  a smoothed  curve  has 
been  drawn  locally.  In  the  surface  conductance  expression,  the 
fin  effect  is  enhanced  by  the  lowered  value  of  0^  as  the  fin 
length  increases  (see  Table  4).  However,  the  temperature  gradients 
along  the  very  long  fin  decrease  in  relative  magnitude.  The  advan- 
tage of  high  flux  from  the  fin  (proportional  to  the  gradient)  has 
reached  a point  of  diminishing  return. 

Similar  calculations  can  be  done  with  fins  in  rectangular 
or  elliptical  ducts  with  anticipated  modifications  from  open  duct 
heat  transfer  values. 

4.  Summary 

The  quantity  2A/bP  = (De/2b)  is  a measure  of  relative  sur- 
face area  that  appears  in  the  definitions  of  Nu  and  Nu'.  It  is 
instructive  to  plot  Nu  and  Nu'/2  (inversely  proportional  to  the 
thermal  barrier  length)  against  the  relative  surface  area  parame- 
ter to  view  the  correlation  of  average  heat  transfer  within  the 
various  ducts;  this  has  been  done  in  Fig.  14.  The  open  ellipsoidal 
curve  essentially  intersects  the  open  rectangular  one,  physically 
reasonable  because  De  approaches  a limit  for  large  eccentricity 
(not  as  in  Ref.  3,  where  a/b  -*•  ® leads  to  a "physically  unrealistic 
limit").  For  open  ducts  of  small  aspect  ratio  a/b,  the  ellipse 
curve  seems  to  be  an  upper  bound,  and  the  rectangle  curve  a lower 
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/2A\  2 2A  De 

\ bP  / bP  2b 


Fig.  14  Correlation  of  Average  Heat-Transfer  Parameters  for  Uniform 

Peripheral  Temperature  in  terms  of  a Relative  Surface  Parameter, 
2A/bP  3 De/2b 
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bound  on  the  ordinary  Nusselt  number  that  can  be  obtained  for  any 
intermediate  shape  (rather  similar  to  conclusions  about  centerline 
velocity  for  the  velocity  field  under  constant  pressure  gradient, 

Ref.  4).  Intermediate  duct  shapes  for  large  a/b  where  De  is 
larger  than  the  elliptic  limit  probably  lie  essentially  on  the  rec- 
tangular curve.  The  apparently  advantageous  high  Nu  with  large  a/b 
is  mostly  due  to  the  Dg  factor  in  its  definition.  The  average  heat 
transfer  based  on  true  geometric  diameter,  Nu',  is  reasonably  con- 
stant for  all  the  open  ducts. 

The  variation  of  Nu  with  a/b  has  been  calculated  elsewhere 
for  open  rectangular  ducts  for  other  homogeneous  thermal  boundary 
conditions  (Refs.  1 and  2).  The  results  have  also  been  plotted  in 
the  coordinates  of  Fig.  14  producing  lower  Nu  curves,  the  one  with 
absolutely  uniform  wall  temperature  having  a similar  trend  to  De. 

The  extended  surface  ducts,  all  derivative  from  the  square 
duct,  are  represented  by  triangles  in  Fig.  14.  Most  of  them,  al- 
though they  have  (2A/bP)  < 1,  have  worse  Nu  than  the  open  square 
duct,  the  data  points  following  the  downward  trend  of  Nu  with  lower 
(2A/bP).  The  average  thermal  barrier  length  is  also  in  the  range 
of  the  open  ducts.  Only  one  point  (starred)  is  on  the  leftward  ex- 
trapolation of  the  upward  trending  rectangular  (Nu'/2)  curve.  Only 
this  finned  duct  is  an  improvement  on  the  open  duct.  This  indicates 
that  merely  decreasing  De  within  an  overall  defined  area  (2a  * 2b) 
is  not  enough  to  ensure  increased  heat  transfer.  An  engineer  must 
be  clever  enough  to  design  the  fin  that  really  is  helpful  and  not 
just  produce  a Nusselt  number  that  follows  the  average  trend  with  De. 


C.  PERIPHERAL  TEMPERATURE  DISTRIBUTIONS 

The  numerical  results  show  that  the  heat-transfer  parame- 
ters are  quite  different  in  a given  duct  with  uniform  peripheral 
temperature  and  in  the  same  duct  with  a distribution  of  temperature 
along  the  perimeter,  even  though  the  average  wall  temperature  is 
the  same. 

The  surface  conductivity  is  proportional  to  the  wall  tem- 
perature gradient.  The  considerable  difference  in  wall  gradient 
when  there  are  wall  temperature  distributions  has  been  shown  by 
the  numerical  calculations.  No  longer  does  the  flux  necessarily 
go  to  zero  when  the  velocity  shear  goes  to  zero.  The  flux  magni- 
tude and  surface  conductivity  are  inversely  proportional  to  6^ 
which  varies  with  the  wall  distribution.  The  variations  of  aver- 
age heat-transfer  parameters  Nu  and  Nu'/2  are  related  only  to  the 
changes  in  9^*,  because  the  other  factors  are  constant.  Because 
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the  distributions  on  the  alternate  square  duct  walls  were  similar, 
proportional  to  e,  the  variation  of  9^  with  e is  linear: 


9b*  ■ <Vo  + eb*(t) 

• -3  J'  u*[80*  + 6*(e)]dA 

Au  A 

1 c * * * 

- J u (%  + E e )dA  . 

Au  A 

60*,  0^*  are  found  in  the  appendix.  Results  at  two  values  of  e 
are  sufficient  to  evaluate  6^  for  all  e. 

The  case  is  similar  for  the  temperature  variation  on  the 
fins  in  a square  duct.  Here  several  values  of  8,  the  value  of  the 
normalized  temperature  at  the  tip  corner,  were  assigned,  and  the 
variation  of  0^*  with  8 is  shown  to  be  linear  in  Fig.  15.  There 
is  a particular  value  of  8 («  0.08)  when  0^  is  zero.  Since  this 
occurs  in  the  denominator  for  Nu  and  Nu'/2,  it  appears  the  aver- 
age heat  transfer  is  infinite  for  this  particular  fin  temperature 
distribution  (and  the  thermal  barrier  length  zero).  A similar  par- 
ticular value  of  e produces  the  same  situation  in  the  open  square 
duct.  Further  changes  in  8 or  e produce  0^  's  of  opposite 
sign.  Physically  this  does  not  make  sense  if  §t/9x  is  a fixed 
(measurable)  quantity;  there  is  either  average  heat  input  to  the 
fluid  per  unit  length  of  duct  or  heat  outflow.  (A  similar  trouble 
was  noted  with  convective  problems  when  the  fluid  contained  heat 
sources  (Ref.  3);  the  situation  was  not  resolved.) 

The  difficulty  seems  to  be  in  the  (arbitrary)  selection 
of  a reference  temperature  (see  p.  4).  Another  convenient  refer- 
ence temperature  might  be  the  centerline  temperature.  Its  varia- 
tion with  6 is  also  indicated  on  Fig.  15.  The  magnitude  is  very 
nearly  the  same  as  0^  • It,  too,  goes  through  zero,  but  at  a dif- 
ferent value  of  8.  The  only  reference  temperature  that  makes  com- 
plete physical  sense  is  the  simple  mean,  0*.  When  that  goes 
through  zero,  the  average  heat  flow  per  unit  length  of  duct  does 
change  sign.  At  that  point  the  perimeter  average  of  the  local 
heat  flux,  as  varied  as  it  may  be  along  the  walls,  is  zero. 

However,  the  practical  interest  is  not  near  zero  heat  in- 
put, but  large  values  of  it.  It  does  not  appear  to  be  a gross 


Bulk  temperature  6 5 
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error  to  use  conventional  0^*  to  compare  duct  configurations  with 
various  temperature  distributions  unless  it  approaches  zero.  It 
certainly  seems  possible  to  increase  the  average  heat  transfer  by 
introducing  wall  temperature  distributions. 

The  actual  temperature  distribution  on  the  duct  walls  is 
due  to  the  thermal  field  in  the  solid (s)  surrounding  the  flowing 
fluid  (Ref.  13).  The  mathematical  problem  should  really  consider 
the  coupled  system,  a conductive  solution  in  the  solid  material(s), 
and  a convective  transport  field  in  the  duct. 


Ref.  13.  E.  M.  Sparrow  and  L.  Lee,  "Effects  of  Fin  Base- 
Temperature  Depression  in  a Multifin  Array,"  Trans.  ASME,  Vol.  97C, 
1975,  pp.  463-465. 
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V.  CONCLUSION 


Enough  examples  of  the  numerical  method  have  been  shown  to 
demonstrate  its  utility  and  to  Judge  its  worth  for  these  particu- 
lar heat-transfer  problems.  The  Cartesian  direct  finite-difference 
convective  solutions  compare  very  favorably  with  accurate  fully 
developed  solutions  obtained  by  other  means.  The  method  is  adapt- 
able for  arbitrary  duct  cross-section  shapes  and  allows  study  of 
various  temperature  conditions.  Normalized  conductance  quantities 
that  have  been  considered  for  heat  transfer  are  shown  to  be  de- 
pendent on  duct  shape  and  perimeter  temperature  distributions  in 
rather  complex  ways,  so  it  is  difficult  to  make  sharp  state- 
ments about  general  trends.  Specific  information  on  wall  tempera- 
ture distribution  is  needed  to  answer  specific  questions  relative 
to  average  or  local  heat  conductance.  Qualitative  information  on 
the  temperature  fields  in  irregular  ducts  can  be  drawn  from  the 
velocity  field  solution,  if  the  wall  conditions  are  similar 
(0W  ~ uw*  * 0).  Close  numerical  estimates  can  be  obtained  from 

the  general  Nu,  Nu'/2  correlations  in  Fig.  14.  But  nonhomogeneous 
peripheral  temperature  changes  the  picture  immensely.  The  examples 
given  point  to  the  inadequacy  of  the  conventional  bulk  (mixed  mean) 
temperature  as  a completely  suitable  reference  temperature. 


The  method  can  be  used  in  conjunction  with  diffusive  field 
calculations  in  solids  surrounding  the  duct,  or  suitable  tempera- 
ture measurements,  to  provide  the  complete  fully  developed  convec- 
tive thermal  field  in  arbitrary  ducts.  The  technique  here  applied 
to  more  or  less  rectangular  shapes  can  easily  be  extended  to  fluted 
annular  ducts,  where  consideration  of  extended  surface  length, 
thickness,  end  rounding,  and  number  per  circumference  can  all  be 
handled  expeditiously  to  optimize  the  convective  field.  Equations 
(la)  and  (6a)  can  be  transformed  to  circular  cylindrical  coordi- 
nates — (n  » in  r * in  (x^/b^  + y2/b2)b),  (<J> ) — where  b is  a 


reference  radius,  <p  is  azimuthal  angle,  and 


-/» L + iL\. 

\ 8n2  3$2/ 


It  should  be  obvious  that  annular  numerical  finite-difference  pro- 
grams would  not  be  very  much  different  from  the  Cartesian  programs 
illustrated  here. 
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Appe  idix 

ANALYTIC  SOLUTIONS  FOR  FULLY  DEVELOPED  TEMPERATURE  FTKLDS  TN  A 
SQUARE  DUCT  WITH  CONSTANT  PEAT  INPUT  PER  UNIT  CHANNEL  LENGTH 


Although  the  details  of  analytic  solution  are  only  given 
here  for  a square  duct,  similar  analyses  can  al90  be  carried  out 
for  arbitrary  rectangles.  Since  the  primary  purpose  is  to  check 
the  numerical  solutions  for  accuracy,  a square  is  sufficient  to 
test  the  computer  programs. 

Normalizing  by  the  central  half-height  of  the  duct,  b, 
and  with  a the  half-width  of  the  duct  and  a » b (square),  the 
velocity  solution  can  be  written  (after  Ref.  A-l)  as 


d£ 
u dz 


k[2  - (X2  + Y2)]  + 3 1 ; 

n*0 


cosh  P X 
n 

cosh  P 


cosh  P Y 


cos  P Y + — 

n cosh  P 


cos  P X 
n 


(A-l) 

2n  + 2.  x v 

where  Pn  * — if , — * X,  and  £ “ Y.  (Note  the  symmetry  in  X, 

Y which  follows  the  physical  expectation  in  a square  duct.)  It  is 
convenient  to  separate  u*  into  two  parts,  the  simple  polynomial 
plus  the  infinite  series: 


u*  - Ul*(X,Y)  + u2*(X,Y)  , 


where  » -1  and  V^U2*  ■ 0,  which  can  be  easily  verified. 

The  solution  to  the  temperature  field  9*  - t*  - t^,*  must 
satisfy  the  normalized  equation 

2 * * 

V 9 - u . (A-2) 


The  solution  to  this  equation  can  be  separated  into  three  parts: 


Ref.  A-l.  V.  O’Brien,  "Pulsatile  fully-developed  flow  in  rec- 
tangular  channels,"  J.  Franklin  Inst..  Vol.  300,  1975,  pp.  225-230. 
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e*  - e^cx.Y)  + e2*(x,Y)  + 03*(x,y) 


(A-3) 


where  V2©}*  - ux*,  V202*  - u2*,  and  V203*  - 0. 

The  first  part  0^*  is  a simple  polynomial: 

* 2 2 2 2 2 
®1  " bl(X + Y ) “ Cj(X  + Y) 

V2©!*  - 4b3  - 16^  (X2  + Y2)  ■ Sj  - k(X2  + Y2) 


(A-3a) 


bx  - 1/8,  c - 1/64. 

Noting  that  V^(Y  sinh  n Y cos  n X)  - 2n  cosh  n X cos  n X,  we  take 


• E °n  ( 


Y sinh  P Y X sinh  P X 

■ — ■ cos  P X + — n cos  P 

cosh  P n cosh  P r 

n n 


* __  / cosh  P \ 

. - V 2P  a ( r-^- 

l n n \ cosh  P 

\ n 


cosh  P X 

cos  P X + -n-  cos  P Y 

n cosh  P n 

n 


(A- 3b) 


) " U2*  ‘ 


so  that 


" 2 P 4 * 

n 


For  0^  we  take  a series  similar  to  u2*: 


* __  / cosh  P H 

3 " 2 6n  (cosh  P 

' n 


cos  P X + — 

n cosh 


P X \ 

p®" cos  v) 

n / 


(A-3c) 


The  coefficients  0n  will  be  chosen  to  satisfy  the  assigned  boundary 
conditions. 
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1.  UNIFORM  PERIPHERAL  TEMPERATURE  tw* 

The  boundary  condition  is  t*  - tw*  * 8*  * 0 at  the  edges 
of  the  square  duct  |Y|  * 1 and  ]x|  - 1. 

(i)  At  Y - ± 1: 


e (x.u)  - 6j  (x,*i>  + e2  (x,±i)  + <x.ti)  - o (All  - l « x < + i) 


n«0 


n*0  n 


n 2 P 
n*0  n 


anh  P cos  P X 
n n 


+ > 6 cos  P X , 

LmJ  n n 

n*0 


(A-4) 


where  use  was  made  of  the  following  identities 


1*2 


£ 

n*0 


(-1) 

P cos  PnX  . 


(A-5) 


1 ~ X2  - 4 ^ cos  PnX  , 


n*0  n 


(A-6) 


2 4 

5 - 6X  + X - 48 


53  cos  P X 

A— ' „ 5 n 

n*0  n 


(A-7) 


tThese  are  taken  from  Jolley  (Ref.  A-2):  Eq"!  (A-5)  from  0506, 
p.  96;  Eq.  (A-6)  from  #528,  p.  100;  and  Eq.  (A-7)  derived  from 
#345,  p.  64. 

Ref.  A-2.  L.  B.  W.  Jollev,  Summation  of  Series.  Dover.  New 
York,  1961. 


- 51  - 


THt  JOHNS  HOPKINS  UNIVt«SlT? 

APPLIED  PHYSICS  LABORATORY 

LAUHfl  MAO 


From  Eq.  (A-4)  we  see  that 


B 

n 


3 (-l)1**1  3 (-l)n 

8 Pn  4 P5 


lidi! 
2 pa 

n 


tanh  P 

n 


(A-8) 


(ii)  At  X - ±1,  the  same  coefficients  also  make  0 “0  because 

of  the  symmetry  in  X and  Y. 

For  easy  comparison  to  the  numerical  results,  we  note  the 
center  value 


6*(0,0)  - 03*(O,O) 


26 


cosh  P 


- -0.06500 


(A- 9) 


n-0 


Another  compact  form  of  analytic  solution  for  a square 
duct  with  the  same  boundary  condition  is  known.  Marco  and  Han 
(Ref.  12)  give  a formula  Involving  double  infinite  summations. 


00 

OD 

, / 2n  + 1 , / 

2m  + 

1 

16  V 

sin  ^ 2 « X | sin  ( 

i 2 

— it  Y I 

* 2^ 

2Li 

\2 

/ 2n  ♦ 1 \ 2 2 

0 

1 

c 

m”0 

(2n  + 1)  (2m  + 1)  \ 

) + 

\ 2 ) 

(A-10) 

where  (X’ ,Y’)  are  measured  from  the  corner  of  the  square  (0  £ X’  t 2, 
0 K Y’  < 2)  by  analogy  to  an  existing  solution  for  deflections  of  a 
simply  supported  thin  plate  under  uniform  lateral  load.  The  velocity 
field  u*  is  also  represented  by  a doubly  infinite  summation. 

The  Nusselt  number  as  a function  of  a/b  (Ref.  12)  is  ex- 
pressed concisely  in  terms  of  double  summations  (A-ll)  and  has 
been  evaluated  by  a short  program  in  APL  language.  (The  author 
is  Indebted  to  Mr.  Stanley  Favin  for  the  program  and  comments.) 
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Nu 


64 


( * 


E 

N-1,3,5,. 


, a,  2 2 

(i  + -)  IT 


mVJm2  + N2  (^) 


) 


1/2  Nu' 


1 


4 (f) 


E 

M 


Nu 


E 

N 


M2N2  M2  + N2(b  ) 


a\2|  3 


(A-ll) 


The  evaluation  of  Nu  (a/b)  with  (K  + 1)  terms  was  easily 
programmed  on  an  IBM  computer  using  the  APL  language.  The  program 
Is  listed  below  along  with  comments. 


HNUCALC 

m 

r 2 1 M*-N+  (l+2*tJO*2 

[3]  nro«-i 

f 4 ] 064*(  (1+j4tB)*01  )*2 
f 5]  MN<-M».*N 


r 6 1 NNA-M* .+N*(A*B)*7 

r 7 1 TOP+( t/ , tMN*MNA '*2 

m BOT*-*-  / , *MN*MNA  • 3 
r 9 1 NU+-C*  TOPiBOT 
riOl  'NU  = ' iNU; ' NU~ / 2 

7 


THE  NAVE  OF  THE  PROGRAM. 

SET  THE  OHIGIN  TO  ZERO. 

GENERATE  TWO  EQUAL  VECTORS  WHOSE  ELEMENTS 

ARE  THE  SQUARED  SEQUENCE  1,3,5,  .(1+2*  ) 
RESET  THE  ORIGIN  TO  1 . 

CALCULATE  THE  CONSTANT  FACTOR: 

C=64/[(  UA/B)**2xPI**2] 

GENERATE  THE  OUTER  PRODUCT  OF  M AND  N. 

MN  WILL  BE  A K*K  MATRIX  OF  NUMBERS  REPRESENTING 
M**2*N**7  INSIDE  THE  DOUBLE  SUMS. 

MNA  WILL  BE  A K*X  MATRIX  OF  NUMBERS  REPRESENTING 
M**2  + N**2*(A /B  ) * * 2 INSIDE  THE  DOUBLE  SUMS. 

TOP  IS  THE  DOUBLE  SUM  SQUARED  IN  THE  NUMERATOR 
OF  THE  EQUATION. 

BOT  IS  THE  DOUBLE  SUM  IN  THE  DENOMINATOR . 

NU 

- ' -,NUx(l+A*B)mxA*B  PRINT  THE  ANSWERS. 

THE  END. 
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2.  BOUNDARY  TEMPERATURE  DISTRIBUTION 

The  differential  equation  is  still  Eq.  (A-2)  but  the  solu- 
tion sought  must  satisfy  6 -eg  at  the  boundary,  where  ?g  - 0. 
Let  the  solution  which  satisfies  0*  - 0 on  the  boundary  (section 
(1)  above)  be  denoted  0q*.  Then  we  seek  0*  - 0q*  h n(e)  where 


V2  n(e)  - 0 (A-12) 


with  B.C.  fig  ■ 0g  . The  symbol  e is  used  because  the  assumed 
symmetric  boundary  distributions  will  be  characterized  by  the  num- 
ber e,  and  the  resulting  solution  will  be  directly  proportional 
to  e. 


At  walls 

|y|  - 1,  let 

0B*(X)  - e(3X2  - 

l); 

(A-13a) 

at  walls 

|X|  - 1,  let 

0B*(Y)  - e(3Y2  - 

1). 

(A-13b) 

The  expression  in  (A-13a)  is  a simple  quadratic  in  X which  can 
be  expressed  as 


0 


*(X)  - (-3e) 


B 

so  if 


n-0  n 


cos  P X + 2e  ~ cos  P X , 

n P n 


0 


i*(t>  ■ E V'>  ( 


cosh  P Y 
n 

cosh  P 


cos  P X + 


n-0 


cosh  P X 
n 


n 


n cosh  P 


cos  P 


n 


"Y)-n 


(e)  , 


E 

n-0 


whence 


cos  P X I 

n I „ 3 P n 

r n 

n . 


B (e)  - 0 , 


Bn(e) 


4c 

P 


(-Dn  + ^4 


n+1  ] 
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By  symmetry  in  X and  Y,  the  boundary  condition  of  Eq.  (A-13b) 
is  also  satisfied. 


* * _ , u; 

The  central  value  0 (0,0)  ■ 0.  (0,0)  + — 2 

0 t—j  cosh  P 

n-0  n 

The  equations  are  linear;  the  additional  solutions  n are  directly 
proportional  to  the  value  of  e.  All  derivative  n quantities  are 
likewise  proportional  to  e;  these  include  not  only  0*(O, 0)  but  bulk 
temperature  0b  , temperature  gradients  (30*/3N)  on  the  wall,  etc^ 
Thus  if  e < 0,  the  central  temperature  can  change  sign  (even  though 
the  average  heat  flux  is  still  positive),  and  even  0b*  can  become 
zero. 

Less  symmetric  assigned  boundary  conditions  on  the  square 
boundary  including  polynomials  of  any  degree  can  be  handled  in 
similar  fashion. 
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1 

LIST  OF  SYMBOLS 

A 

2 

Cross-sectional  area  of  duct  (m  ) 

, 

a 

Horizontal  semidiameter  of  duct  section  (m) 

b 

Vertical  semidiameter  of  duct  section  (m) ; characteristic 
length  (Fig.  1) 

c 

P 

Specific  heat  of  fluid  (W*s/g,°C) 

d 

Depth  of  fin  (m) 

D 

e 

Hydraulic  diameter  (m) 

h 

Finite-difference  mesh 

h 

2 

Thermal  conductance  (W/m  , °C) 

k 

Thermal  conductivity  of  fluid  (W/m,  °C) 

c 

N 

Normal  direction  to  wall  surface 

Nu 

Conventionally  defined  Nusselt  number  (p.  33) 

IJ 

Nil’ 

Unconventionally  defined  Nusselt  number  (p.  34) 

fr 

P 

Duct  section  perimeter  (m) 

P 

2 

Pressure  (dyn/m  ) 

[ 

Q 

2 

Volume  heat  addition  per  unit  channel  length  (W/m  ) 

s’ 

Diagonal  stretched  coordinate  (p.  22) 

j I] 

t 

Temperature  (°C) 

ii 

* 

t 

Normalized  temperature  (Eq.  (6b)) 

i; 

u 

z-ward  velocity  component  (m/s) 

li 

* 

u 

Normalized  velocity  component  (Eq.  (lb)) 

t 

x.y.x 

Cartesian  coordinates  (m) 

l 
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X,Y  Normalized  Cartesian  coordinates  (xb  \ yb  ^ , respectively) 

a Thermal  diffusivity  of  fluid  (m/s) 

6 Tip  fin  corner  temperature  (normalized) 

c Parameter  characterizing  peripheral  wall  temperature 

(Eq.  (A-13)) 

n Convective  temperature  field  due  to  peripheral  wall  tem- 

perature distribution  (appendix) 

* * * 

0 Normalized  temperature  difference  (t  - t ) 

w 

V Viscosity  coefficient  (g/m  s) 

3 

p Fluid  density  (g/m  ) 

Superscripts 

( ) Sectional  average 

* Normalized  quantity 

* Modified  coordinate  (or  parameter) 

Subscripts 

b Bulk  (or  mixed  mean)  value 

B Boundary  value  (appendix) 

c Central  value 

w Wall  value 

ref  Reference  temperature 

loc  Local  (surface)  value 

m Summation  index  (appendix) 

n Summation  index  (appendix) 
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